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1.0 INTRODUCTION 


This document is the final report covering the Satellite Borne Gravity Gradiometer Study con- 
ducted by Bell Aerospace under Contract NAS 5-20910 to Goddard Space Flight Center. 

A number of earth orbital gravity measuring techniques have been studied and conducted in the 
past. Satellite-to-satellite doppler velocity measurements and groimd tracking of satellite orbital perturb- 
ations have successfully yielded the long wave length gravity anomalies of the earth. These techniques 
give useful gravity anomaly data below the 20th harmonic degree but little information at the higher ones 
because the orbital perturbations are related to the gravity anomalies through an integration process. 

Orbital gravity gradient measurements on the other hand are related to the gravity anomalies 
through spatial differentiation. Even though the gravity gradients are very much attenuated at orbital 
altitudes, the analysis discussed in Section 2.0 of this report shows that gravity anomaly determination to 
the 100th harmonic degree is possible with an orbital gravity gradiometer in the 0.03 EU (10 sec avg. 
time) (1 Eotvos Unit = 10“®/sec^) performance category. 

Gravity gradiometry is recognized to be a very difficult instrumentation problem because extremely 
small differential acceleration levels have to be measured, 0.1 EU corresponds to an acceleration of 10‘“g 
at two points 1 meter apart. Bell Aerospace is under contract to SAMSO USAF since April of 1974 to 
develop a feasibility model of a gravity gradiometer for airborne applications using four modified versions 
of the proven Model VII accelerometers mounted on a slowly rotating fixture. It is an objective of the 
program to demonstrate gravity gradient measurements to an accuracy of I EU in a ten second time period 
in all orientations with respect to the g vector and under moving base conditions in the laboratory. The 
progress so far has been excellent and gravity gradients are being measured to 1 .07 EU in a vertical rotation 
axis orientation. Equally significant are the outstanding operational characteristics such as fast reaction 
time, low temperature coefficients and high degree of bias stability over long periods of time. The rotating 
accelerometer gravity gradiometer approach and its present status is discussed in Section. 4.0 and is the 
foundation for the orbital gravity gradiometer analyzed in this report. The performance levels achieved at 
this time, in a 1 g enviomment of the earth and under relatively high seismic disturbances, lend the orbital 
gravity gradiometer a high confidence level of success. 

A slowly rotating satellite (0.5 to 0.1 rad/sec spin speed) at an altitude of between 250 and 350 km 
would be an almost ideal vehicle for a gravity gradiometer. The low g and vibration free environment and 
the “smooth rotation” provided by a spinning vehicle will eliminate many of the error mechanisms present 
in the airborne gravity gradiometer. Furthermore the low g environment will permit the use of the Space 
qualified MESA (Miniature Electrostatic Accelerometer) with the inherent lower noise. The electrostatic 
suspension system of the MESA system provides near zero spring constant along the sensitive axis and 
eliminated one of the major noise producing sources stemming from- the spring suspension of the Model 
VII proof mass. Over 50 MESAs have successfully operated on various-orbital missions including: 

No. of Missions Name No. of MESAs per Mission 


1 

2 

1 

7 

1 

1 

3 

3 


Logacs - Agena 
Cannonball Satellite 
Spades 
Agena 110 
SERT II Satellite 
S-73-5 SatelUte 
Atmospheric Explorers 
DMA -501 Experiment 


1 

3 

1 

1 

1 

1 

3 

1 three axis MESA 


1 



The only modification contemplated to the MESA for the orbital gravity gradiometer application 
is hard evacuation of the instrument to eliminate gas damping. Hard evacuation was successfully carried 
out on the Model VII accelerometers and has been maintained for over a year on five instruments. 

This report is organized into sections which develop a logical case for an orbital gravity gradio- 
meter. 


The analytical section covers a discussion of: 

(1) The statistical nature of the gravity anomalies expected to be encomitered at an orbital 
altitude of 300 km and their relationship to gravity gradient anomalies. 

(2) An estimation of accuracy with which the gravity anomalies will be measured as a function of 
gravity gradiometer performance and harmonic degree of the earth. 

(3) An estimation of the accuracy with which the gravity anomalies will be me^ured as a function 
of satellite tracking accuracy and harmonic degree of the earth. 

A descriptive design section of the orbital gravity gradiometer which covers; 

(1) The description of the rotating accelerometer orbital gravity gradiometer outlining the tech- 
niques of this concept which permits high accuracy gravity gradient measurements with 
proven and qualified instruments. 

(2) A description of the Miniature Electrostatic Accelerometer 

(3) Discussion of MESA instrument noise and recommended modifications to reduce it to 
required level 

(4) Discussion of orbital vehicle induced noise and outlining orbital vehicle critical design para- 
meters for the gravity gradiometer mission 

(5) A status report of the airborne rotating accelerometer gravity gradiometer development 
program, which forms the basis for the orbital gravity gradiometer 

(6) Short section summarizing the conclusions and recommendations. 
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2.0 DISCUSSION AND ANALYSIS OF ORBITAL GRAVITY GRADIENT MEASUREMENTS 


Satellite gradiometry offers the potential of determining the global gravity field of the earth with 
greater detail than alternative satellite schemes. The investigation reported here explores the relative merits 
of a gradiometer system and a doppler velocity system. The two systems are complementary in that the 
doppler velocity system is superior for defining the lower harmonic degree variations (less than 30th har- 
monic degree). The gradiometer system on the other hand gives superior performance for the higher har- 
monic degrees. The investigation proceeds through four major problem areas: 

(1) Determination of a statistical characterization of the sea level gravity anomaly of the earth. 

(2) Determination of statistical characterizations of various elements of the gravity field and satellite 
perturbations at satellite altitude. This statistical characterization is in the form of Power 
Spectral Density (PSD) curves for the anomalous gravity gradient, gravity anomaly, and finally 
the radial velocity perturbations of the satellite. The last of these provides a basis for evalua- 
ting the doppler velocity system. The procedure followed here is somewhat in contrast with 
the usual one of using Kaulas rule of thumb to characterize the -various anomalies and velocity 
perturbations at satellite altitude. As will be shown, however, the results are in good agreement 
particularly at the higher harmonic degrees. 

(3) Resolution of instrumentation and associated problems such as satellite attitude, altitude, and 
position control, optimal orbit determination data handling, etc. 

(4) Finally the optimal estimation problem which involves a parameterization of the earth’s 
gravity field and utilization of satellite measurements to establish parameter values. 

A discussion of each of these general problem areas is presented to outline the approach, the anti- 
cipated problem areas and expected results. 

2.1 DETERMINATION OF A STATISTICAL CHARACTERIZATION FOR THE GRAVITY FIELD 

AT THE EARTH’S SURFACE 


A statistical characterization of various gravity related signals (gravity gradient and gravity anomaly, 
satellite velocity perturbations etc) at satellite attitude is central to an assessment of various satellite based 
geodetic systems. The usual approach is to use Kaulas’ rule of thumb to define these signals. The approach 
taken here is to use a sea level model for the gravity anomaly and to analytically extend this model to 
satellite attitudes. Extensive sea level gravity anomaly survey data provides a basis for this approach. 


A number of investigations (References 1,2,3, 4) have dealt with this problem. The model 
illustrated by the block diagram on Figure 1 is generally consistent with most of the models developed. A 
second, shorter correlation distance (Dj ) has been introduced so that the sea level gravity gradient anomaly 




component Tj^^ ^ ) is 


^x-^z 


bounded and in agreement with RMS levels expected. 
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Figure 1 . Earth Surface Model of Gravity 


where W = U + T 
W - Gravity potential 

U - Normal component of gravity potential 

T - Anomalous component of gravity potential 

- i = X, y, z gravity anomaly components 

Tjj - i, j = X, y, z anomalous components of gravity gradient 

- 20 nm, 2.5 nm corelation distances 

ot- - 40 mgal 

2.2 RELATIONSHIP BETWEEN SURFACE ANOMALIES TO MEASURABLE VALUES AT 
ORBITAL ALTITUDES 

At satellite altitudes, the higher spatial frequencies present in the sea level signals will be greatly 
attenuated. Pizzetti’s extension of Stokes integral provides the analytical basis for extending the sea level 
model to satellite attitude. This integral equation relates gravity potential at satellite altitude to gravity 
anomaly on the earths surface. The following development details how the gradiometer outputs, which 
consist of various second spatial derivatives of gravity potential, are related to the sea level anomaly. 

A rotating orbital gravity gradiometer in a polar orbit, whose sensitive plane is aligned to the orbit 
plane, has two outputs . A difference of in line gradients and a cross gradient as follows: 

E, - 7 )W .X) (= - T,,,,) ( 1 ) 


Ej (r,0,X) ^ T(r,0,X) (= IT ^ 

r 9r90 

where - latitude, longitude 

r - orbital radius of an assumed circular orbit 

* 

T(r,^,X) - anomalous gravity potential at orbital altitude 


( 2 ) 
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Pizzetti’s extension of Stokes integral relates T(r, 0 ,X) to earth surface anomalies; 

T(r,0,X) cos<^'d<^'d\' 

where S is Stokes function given by 

. r 2 / 1-5 cosi// + D 

S(r, 0 ,X, 0 ' ,X' ) = — + 1 - 3D - 1 cos ^5 + 3fin (■ ^ — 

and cosil' = sin^ sin 0 ' + cos^ cos ^ ' cos (X-X ' ) 

t = R/r 

D = 1 - 2t cosil' + 1^ 

R - mean radius of earth 

Substituting 3 in 1 and 2 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


^ R {d^S. a^s 
~ 4n ^ I 3r2 ' r^ 3^2 


5gCOS0' d^'dX 


R r 2 3^ S 1 


( 7 ) 

( 8 ) 


The bracketed terms in (7) and ( 8 ) can be viewed as weighting functions which relate gradiometer 
outputs to anomaly at ^ ^ , X V Normalized curves for the bracketed terms in equations (7) and (8) are 
plotted in Figures 2 and 3, respectively. The solid curves in Figure 2 show the weighting function for three 
satellite altitudes as a function of along track displacement in degrees. The dotted curves show variation 
with cross track displacement in degrees. The sharper more localized character of the weighting function 
for low satellite altitudes translates into higher resolution. A spinning orbital gravity gradiometer yields 
lower resolution cross track than along track. The weighting function shown in Figure 3 
23^8 

(rgra^) for the gradiometer output £3 has odd symmetry with respect to along track displacement. E 2 

must approach zero mean over long orbital distances. During a satellite pass over a discrete localized sur- 
face anomaly it is evident from the weighting functions that the peak to peak outputs of E 2 and Ei are 
the same. Therefore as much information about the earth’s gravity field is contained in E 2 as in the signal 
Ej . Power spectral density plots will show this more explicitly. This result is in contrast to Reeds [7] re- 
sults which show much lower signal levels for T ^2 Ty^. 

I — 2 2 3T 

Weighting functions for gradiometer amplitude output E, (=\/Ei + E 2 ■ ) and for Ty (=-ry are 

^ or 

shown in Figure 4. Both of these are seen to be less localized than Ej and therefore we should expect 
lower resolution of the gravity field based on these measurements. Since radial velocity perturbation is 
dynamically the integral of we conclude that the weighting function for 5 V is even less localized than 
that shown for T^. 

The highly localized nature of the weighting functions for Ej and -Ej permits an approximation 
in equations (7) and ( 8 ) causing them to take on the form of two dimensional convolution integrals. Con- 
sequently Fourier transform techniques can be used to characterize these relationships. In particular, 
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Figure 2. Normalized 
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Figure 4. Normalized Weighting Functions for Ej, E 3 , 

~ for 300 km Altitude 


'frequency domain characterization of the signals Ej and E 2 are possible. To complete the picture all that 
is required is to relate velocity perturbations of the satellite to these signals. The following matrix differen- 
tial equation relating anomalous gravity at satellite altitude provides the connection. 


T T 
^ XX xy 



T T 
^xz .^yz 



T V 
^xx '^x 


T V 
^xy ''x 


T V 
^xz ''x. 


(9) 


In writing equation (9) a circular orbit has been assumed with the x axis directed along track. The dynamics 
of satellite velocity perturbation is quite complex, however, for the higher harmonic degrees and for error 
analysis purposes the simple dynamics depicted in Figure 5 is an adequate model of the radial dynamics. 



Surface 


Figure 5. Radial Dynamics Model 
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This model.relates 6g (earth surface anomalies), the gradiometer output and radial velocity perturba- 
tion 5 V^. A similar model relating Ei and along track velocity perturbation is not directly possible since 
El is a measure of in linear combination with Tyy . 


Using the earth surface model of gravity anomaly and the Fourier transform techniques, the gravity 
gradient power spectral densities for the signals Ej and E2 have been calculated. The power spectral den- 
sities for a 300 km orbit are shown in Figure 6. Glaser and Sherry [6] using Kaulus’ rule of thumb for 
T22 (equ^t.ion 10) have approximated the signal Ei (= T^^ - by 

[El (n)P« [2T22(n)]2 ^\— k 10®J f- j (n+ 1) (2n + 3)^ (2N+ 1) 

( 10 ) 

where n - harmonic degree. 


Tlris’ curve is also plotted in Figure 6 along with Reeds [7] estimate of the gradiometer amplitude output 
(E3 ) which is also based on Kaulas’ rule of thumb. 

E," - E,’ + e " - ClI) 


All of these estimates for gradient signals at orbital altitude are in fairly close agreement particularly at the 
higher harmonic degrees. 


Figures 7, 8 and 9 summarize the PSD’s of various signals at three circular orbital altitudes. Table 1 
summarizes a set of pertinent values for the three altitudes. Table 2 summarizes the RMS statistics for 
the signals E2 and T2 . 


Assuming that the gradiometer measurements and the radial velocity doppler measurements are 
limited by additive white noise, conclusions can be drawn about the required white noise levels and the re- 
lative merits of the two instrumentation schemes. Note that white noise of the gradiometer, which has a 

constant PSD in the Ej^ , E2^ plots (Figure 7) varies asf— on the plot (Figure 8) and as(— M' 

\HD/ ^ \HD/ 


on 


the 6 plot (Figure 9). By the same token white doppler velocity error (constant PSD on S V’ plot 
Figure 9) varies as (HD)^ in the domain (Figure 8) and as (HD)^ in the Ej^ , £2^ domain (Figure 7). 

Three levels of gradiometer white noise and three levels of doppler velocity white noise are superimposed 
on the PSD plots of Figure 7. A measure of the resolving power of a particular instrumentation scheme is . 
the harmonic degree at which the instrument white noise level is equal to the signal PSD. Table 3 sum- 
marizes the results for the gradiometer while Table 4 does so for the doppler velocity instrumentation. 

This comparison is based on using the gradiometer output E2 and the radial velocity perturbation 6 V^. The 

additional value derived by utilizing Ej has not been assessed nor has that of instrumenting along track and 
across track velocity perturbations. Figure 7 shows that for any given white noise levels for the gradiometer 
and the doppler velocity system there exists a crossover harmonic degree below which the doppler data 
is of more value than the gradiometer and above which the reverse holds true. Since resolution is a function 
of the highest harmonic degree of signal power that can be measured, we see that the gradiometer in- 
herently is capable of higher resolution providing the instrument can be built with a sufficiently low white 
noise level. 
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Figure 9. Power Spectral Density Curves for 






TABLE 1 

SATELLITE AND GRADIENT VALUES 


Orbit Altitude 
(km) 

GM / 1 \ 
r3 

0 

u 

(r/s) 

1 GM \ 
^ ) 

Orbit Period 
(Min.) 

V (m/s) 

0 

1.536E-6 

0.00124 

84.49 

7910 

200 

1.400E-6 

0.00118 

88.49 

7760 

300 

1.338E-6 

0.00116 

90.52 

7750 

500 

1.225E-6 

0.001 1 1 

94.62 

7630 


GM = 0.398603E15 m^sec^ 
R = 6.37816E6 m 


TABLE 2 

RMS VALUES FOR E2 AND T^ 


Orbit Altitude 
(m) 

e2 

Tz 


ff(E) 

a^(mgal^) 

ff (mgal) 

200 

0.0998 

0.315 

66.6 

8.2 

300 

0.0360 

0.190 

48.2 

6.9 

500 

0.00812 

0.090 

.26.9 

5.2 


TABLE 3 

GRADIOMETER INSTRUMENTATION 


White Noise Level 

Harmonic Degree of Resolution 
at Orbital Altitude 

E^/HD 

EVr/s 

E in 10 sec 

200 km 

300 km 

500 km 

icr^ 

0.847 

0.291 

57 

36 

-- 

1CT‘* 

0.0847 

0.0921 

111 

67 

36 

10“^ 

0.00847 

0.0291 

150 

96 

54 


TABLE 4 

DOPPLER VELOCITY INSTRUMENTATION 


White Noise Level 

Harmonic Degree of Resolution 
at Orbital Altitude 

(M/S)^ 

HD 

(M/S)2 

r/s 

cm/s in 
10 sec 

200 km 

300 km 

500 km 

2.95 X ICr® 

2.5 X 1CT® 

0.05 

48 

38 

28 

2.95x1Cr‘° 

2.5x 1CT’ 

0.0158 

70 

54 

38 

2.95x1CT“ 

2.5 X 1(T® 

0.005 

92 

69 

48 







At the 300 km altitude a 6.9 mgal using a 0.03E in 10 sec gradiometer can be estimated 
to a level of aT^ = 0.66 mgal. Using a 0.005 cm/s in 10 sec doppler the result is - 0.12 mgal. If both 
instruments are used the result is = 0;084 mgal. 

2.3 INSTRUMENTATION PROBLEMS 

The instrumentation requirements associated with determining satellite position over the earth, its 
altitude, attitude, etc. for the doppler velocity system differ from those for the gradiometer system. A 
discussion of some of the salient aspects of these requirements for the two systems follows. 

2.3.1 Doppler Velocity System 

From a practical point of view, a doppler system involves two or more satellites; one of them being in 
a relatively low altitude polar orbit, the second being in a higher orbit. The doppler system is instrumented 
to measure the velocity along the line of sight between the two satellites. Several problems must be re- 
solved. 


(1) The velocity perturbation attributable to the low satellite must be determined. The measured 
signal includes the total velocity difference along the sight line. 

(2) The velocity perturbation of the low satellite attributable to gravity anomalies has to be ex- 
tracted from all sources causing perturbations including solar pressure and orbital drag. 

(3) The orbital parameters of the low satellite are varying continuously. This fact presents a 
problem as to which set of orbital parameters the velocity perturbation is referenced to. This 
argument applies to the high satellite as well, but its ephemeris changes more slowly. 

2,3.2 Gradiometer System 

By using a set of four MESA’s the gradiometer is capable of discriminating against atmospheric 
drag, solar pressure forces as well as magnetic torques which may act on the satellite. As a result, the data 
reduction problem is greatly simplified. The gradiometer mission can be accomplished with a single low 
altitude polar orbit satellite. A fairly uniform set of measurements covering the globe can be collected in 
a matter of days. 


2.3. 2.1 Satellite Altitude Control - If the signal Ej is to be used, an orbit altitude determination require- 
ment must be imposed. The normal gradient (based on a spherical earth) is given by: 


Ua 



-1 0 0 \ 

0 -1 0 ) 

0 0 2 / 
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Variation with, altitude is given by 


aUy -3GM 


at 300 km 


9r 
GM 


= 1338Eand 



= 7.23 X 10-'* E/m 


So the normal component of Ei attributable to radial error is 


— = 3 (^^) = 0.00217 E/m 
dr - r 


If dEi is to be held to 0.03E then dr must be known to 13,8 meters. 

2.3. 2.2 Attitude Control Requirements - The major component of the gradient tensor, the normal com- 
ponent, is that due to the earth considered as a point mass. More generally the normal component may 
represent the gradient due to a low order spherop. The anomalous component used in the data reduction 
is computed by subtracting from the gradiometer outputs the appropriate normal components. The 
satellite coordinate system and the computer coordinate system in which the normal component is com- 
puted must therefore coincide to some accuracy. The required accuracy of satellite attitude control or 
determination is analyzed as follows. 


The local level orbit coordinate system has its X axis directed along the velocity vector Y axis 
(satellite spin axis) perpendicular to orbit plane and Z axis vertically downward. The normal gravity 
gradient in this system is 





0 

-1 

0 


0 ^ ™ CUN) 


( 12 ) 


For the more general case the direction cosine matrix relating the normal coordinate system to the 
satellite coordinate system can be defined in terms of a sequence of angle rotations 


( cosi/i sini/c 0 \ 

-sulci' cose/' 0 J 

0 0 1 / 



0 

COS0 

-sin^ 



cose 

0 

-sine 

0 

1 

0 

sine 

0 

cose 


A small angle deviation from this nominal orientation can be treated as a vector. 


( 13 ) 




1 

8^ 


1 

56 

-S(j> 

= I 

+ [6A] 



/ 1 0 0 \ / 0 5 ^ 

f 0 1 0 ] + ( h5(|' 0 

\ 0 0 1 / \s8 -8(p 



( 14 ) 
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The error in gravity gradient ATjj due .to the attitude uncertainty is given by; 

(I + [A A] ) C-/ [UN] (%^) (I + [5 A] )T - fUN] (C^S) j ( j 5 ) 


ATij 


ATn = 


r-” 

GM 


r3 


Cn^ [UNI (Cn^)"^ [ 6 A]T+ [5 A] %S [UN] (Cn^)T 


+ [5 A] [UN] [5A]T 


(16) 


For the special case where the satellite is controlled to the nominal attitude = I and equation 1 6 
becomes 


AT;; = 


GM 


[UN] [5A]T+ [5 A] [UN] + [5 A] [UN] [ 6 A]'^ 


ij j... 

Evaluating this expression, the result to first order off diagonal and 2nd order diagonal is; 

-5}p^+2SB^ 0 -386 


(17) 


AT;; 


GM 


0 


-360 


-6 0^ + 2 60^ 350 

360 -50^ - 50?/ 


(18) 


The conesponding error in the two gradiometer outputs are therefore; 

GM 


GM 

AEj = — 


-502 -602 -(- 502 + 2502 ) = 


-3502 .§02 ^§02 


^ GM ^ ^ ^ ^ GM ^ 

AE -2 ='— 2 -350 =-—660 

r^ r2 

GM 

Equation 20 implies A Tjj2 “ — 3“^^^ 


(19) 

( 20 ) 


As described elsewhere in this report the signal Ej can be used on a long term average basis to define the 
pitch attitude of the satellite. As a result 60 can be controlled or known to negligibly low level. If the 
error in E; is to be controlled to less than 0.03 E for example, equation (19) requires, 


602 


502 ^ 


0.03 

GM/r2 


— — — - 2.24 X ICT® rad2 
1338 


( 21 ) 


for the 300 km orbit. If we neglect the higher harmonic degree attitude perturbations the satellite axis 
tends to remain stationary in inertial space, hence roll and yaw errors vary at orbit rate- as follows: 


50 = 60 sin (w^t + T) 
60 = So cos(Wgt + 7) 
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where Sq is the total angular deviation of the satellite axis relative to the orbit plane direction, 

(Uq - orbital rate 

7 - defines the direction of the angular deviation 5 q 

Evaluating 21 the result is: 

5,^2 _g ^2 = cos ( 2 Wq 1 ) 

If we require the peak error < 0.03E then Sq. < 0.27 1“. 

With respect to determining position, altitude, attitude, etc. the potential of accomplishing real 
time space navigation merits investigation. A satellite equipped with a gradiometer made up of sensitive 
accelerometers such as the MESA, has the basic elements necessary to implement real time space naviga- 
tion. The low polar orbit satellite will experience significant drag and gravity anomaly forces which can be 
sensed by the MESA’s and the gradiometer. Real time ephemeris updating should be possible using this 
data. The gradiometer data is also useful in establishing a coordinate system by virture of defining the 
local vertical. Real time orbital parameter computation could lead to significant relaxation of ground 
tracking requirements'. 

2.4 DETERMINATION OF THE EARTH’S GRAVITY FIELD BASED ON MEASUREMENTS 

MADE IN ORBIT 

At the outset, it would appear that a spherical harmonic parameterization of the earth’s gravity 
field would be the most natural. Several investigators [5,7], however, have used an alternative parameteri- 
zation which divides the earth’s surface into lattitude, longitude squares. The optimal estimation problem 
involves determination of the mean anomaly for each of these areas. This alternative parameterization is 
made possible by the highly localized nature of the gravity gradient weighting functions. Argentiero and 
Garza-Robles [5] demonstrate near orthogonality of a core of squares imbedded in a larger block of gradio- . 

meter data. The spherical harmonic parameterization by contrast does not permit localized solutions based 
on localized data. 

The long term average output of the gradiometer signal Ej (= 2Tj. must statistically approach zero. 
The output Ej (= W 22 ■ other hand includes a normal component of~^—(=» 4014E at 300 km). 

Consequently Ej or E 3 (Gradient amplitude) measurements to an accuracy of 0.03E require scale factor 

stability of 7.5 x lO”® Scale factor stability to this level may be difficult to achieve. The impli- 

cations of this as it relates to the data reduction problem are pervasive. A geographically localized block of 
data in general spans several days in time. If the signals Ei or E 3 are to be used in the data reduction it will 
be necessary to model and estimate or control long terra variations in scale factor. The net result is that 
cross track resolution will be further degraded over that of the basic instrument. The long term stability 
problems might be controlled by utilizing the high redundancy of data taken over the polar regions. Short 
term signal variations (IQth to lOOth harmonic degree) occurring during a single pass will be much less cor- 
rupted by scaleTactor variations. The signal being a zero mean signal is virtually unaffected by scale 
factor stability. 
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A detailed analysis of the data reduction problem utilizing the signals Ej and E 2 and modeling 
long term scale factor stability should be carried out to resolve these problems. Such an analysis would 
also answer quantitatively the questions of cross track resolution and along track resolution for estimating 
the earth’s gravity field. 
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3.0 ROTATING ACCELEROMETER GRAVITY GRADIOMETER CONCEPT 


The previous section outlined the merit of global gravity mapping using an orbital gravity gradio- 
meter. Gravity anomalies could be mapped to much higher harmonic degrees and in a shorter period of 
time than possible with present techniques. The question therefore is whether an orbital gravity gradio- 
meter is feasible within the present state of the art of instrumentation, at a realistic cost and with sufficiently 
high confidence level of success to warrant a satellite mission for this purpose. 

Bell has been under contract to SAMSO USAF to develop a feasibility model of a rotating acceler- 
ometer gravity gradiometer for airborne application since April of 1974. It is the objective of that pro- 
gram to solve the very difficult instrumentation problem of moving base gravity gradiometry within the 
present state of the art and application of proven and qualified instruments. Sufficient progress has been 
made to demonstrate gravity gradient measurements in the laboratory to 1 EU in a ten second time period 
with excellent operational characteristics. 

3 . i DESCRIPTION OF ROTATING ACCELEROMETER GRAVITY GRADIOMETER CONCEPT 

^Dynamic gravity gradiometry constitutes an extremely difficult instrumentation problem. 0.03 EU 
gravity gradient corresponds to an acceleration difference of 3 x 10"‘^g at two points 1 meter apart which 
has to be measured in the presence of vehicle disturbances. On the earth this will include vehicle acceler- 
ation in addition to the full 1 g gravity vector. In orbit the maximum acceleration levels experienced will 
be under 10'^ g. This much lower acceleration environment in space is the major reason that gravity 
gradient measurements to 0.03 EU are deemed feasible aboard an’ orbital vehicle. 

At first glance measuring acceleration differences to 10‘*^g at two points 1 meter apart appears to 
pose an unsurmoun table instrumentation problem. 

The null bias, the output in the absence of an input, of a good force rebalance accelerometer, such 
as the Bell Model VII, is in the order of 20 micro g with a randomness of a few tenths of microg’s as illus- 
trated on Figure 10. At first glance extracting acceleration inputs 5 order of magnitudes smaller than the 
null bias noise does not appear feasible. A power spectral analysis of a typical null bias run of many hours 
is shown in Figure 1 1 . The power spectral density of the null bias is high at the low frequencies but falls 
off to very low level at frequencies above a revolution per minute. 

Null bias is produced by non-acceleration induced forces on the proof mass of the accelerometer. 
Figiue 12 shows a cut-away section of the Model VII accelerometer which is used with the airborne gravity 
gradiometer, to illustrate the nature of null bias and also to highlight the benefits obtained by using the 
Miniature Electrostatic Accelerometer for the orbital gravity gradiometer. The cylindrical proofmass is 
shown suspended by flexure springs. The proof mass has a tendency to hang back with acceleration along 
the sensitive or input axis. This deflection is sensed by a capacitive pick-off. An electronic loop operating 
from this pick-off signal generates a current which is injected into a coil mounted inside the proof mass. 

This current interacts with a strong magnetic field produced by the conical Alnico V magnets to restore 
the proof mass to the null position. The force rebalance current is the output signal of the accelerometer. 

Should a plate of the capacitance pick off system be slightly displaced or should the zero torque posi- 
tion of the suspension spring change slightly because of a change in temperature or aging characteristics of 
materials, a current would be required to maintain the proof mass in the pick off null position. For these 
reasons all inertial instruments are to some extent thermometers and indicators of aging characteristics of 
materials and hence the slow variation in null bias with time and high bias power spectral densities at very 
low frequencies. 
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Figure 12. Model Vll Accelerometer 


The MESA with, its scalable electrostatic suspension system exhibits much better null bias charac- 
teristics than accelerometers with conventional flexure suspension systems. This is discussed in the next 
section of this report. 

If the gravity gradient information could be shifted to a higher frequency region in the order of a 
revolution per minute or more, detection of the minute gravity gradient' signals would clearly become 
possible as indicated by the power spectrum in Figure 1 1 . This is accomplished by mounting four acceler- 
ometers on a rotating fixture as illustrated schematically by the diagram of Figure 13. Each accelerometer 
will still exhibit its own null bias in the very low frequency domain. 

The input linear acceleration experienced by each accelerometer will be modulated at the rotation 
speed n , and the gravity gradient signal will be modulated at twice rotation frequency 2n • The gravity 
gradient signal has been separated in the frequency domain from the unwanted noise by O and 2J2, and 
the information can be extracted by synchronous demodulation. In addition, unwanted signals are can- 
celled by common mode rejection made possible by proper summation of the four accelerometers. The 
output of opposing pair of accelerometers (aj a 2 ) and (as + a 4 ) are summed to reject the linear acceler- 
ation within the scale factor balance of the two instruments for two reasons. Since the fixture is rotated 
relatively slowly it would require a very long filtering time to separate the large output at 1S2 caused by 
linear accelerations from the very small gravity gradient output at twice rotation speed 2Cl . Furthermore, 
linear acceleration inputs at or near the rotation frequency would produce an error signal at twice rotation 
frequency. 

A spinning orbital vehicle will experience some linear acceleration inputs at twice rotation speed 
because of orbital drag being modulated by imperfections in the geometry of vehicle. A second pair of 
accelerometers is mounted to the fixture in space quadrature and their summed output subtracted from 
tlie first to minimize the effect of angular' modulation at twice spin speed wliich would be fully measured 
by only one pair. An additional benefit of the four accelerometer is the multiplication of the gra'vity 
gradient output by a factor of two. 


The rotating accelerometer gravity gradiometer is inherently a two axis device, measuring the 
difference of in-line gradients on one channel and the cross gradient on the other. For example, a rotating 
accelerometer gradiometer with a horizontal rotation axis moving along the north (x) axis is illustrated 
schematically by the gradiometer in the background of Figure 13, and this will be the orientation of the 
satellite gradiometer in a polar orbit except of course at much higher altitudes. 

3.2 ACCELEROMETER SCALE FACTOR BALANCE LOOPS 

As pointed out, pairs of accelerometers are used with the rotating accelerometer gravity gradiometer 
to reject linear accelerations within the scale factor balance of accelerometers in each pair. If the scale 
factors of the accelerometers in each pair were exactly equal, linear acceleration inputs would be totally 
cancelled because the input axis of the accelerometers in each pair are pointing in opposite directions. Long 
experience with inertial instruments has shown that balances between two accelerometers to 1 0"'^ are 
relatively easy to achieve and maintain. Balances of 10’® or better are not only very difficult to attain but 
cannot be maintained very long because of material and electronic component instabilities. To tolerate a 
reasonable acceleration environment at the spin frequency, scale factor balances of lO"^ are desirable. 

However, the system is continuously informing us as to the state of the scale factor balances of the 
two accelerometers in each pair by an output signal at the rotation frequency in conjunction with a steady 
state acceleration input. On the earth the steady state acceleration is provided by gravity and in orbit by 
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orbital drag. The output of the gradibmeter is demodulated by the sin and cos of the spin speed (ii ) and 
the resultant signal integrated and applied to adjust the scale factor of one accelerometer in each pair until 
the signal at the rotation frequency disappears. 

Block diagram Figure 14 illustrates the scale factor balance loops. Theoretically a very higli degree 
of scale factor balance is achievable because relatively long time constant loops are permissible since scale 
factors change only slowly. With the Model VII accelerometer the scale factor balances are. adjusted by 
injecting a current into coils woimd around the Alnico magnets. The MESA scale factors are adjusted by 
controlling the dc forcer voltage. These scale factor balance loops have been implemented and are oper- 
ating well. 

The second scale factor balance loop adjust the summed scale factor balances of accelerometers 
aj + aj = aj + a 3 to reject angular modulation of the spin speed at twice rotation speed. A considerable 
amount of twice speed angular modulation has to be expected from the rotating mechanism, bearings, 
torque motor, tachometer and slip rings of the airborne gravity gradiometer. The rotating fixture is 
oscillated through a small angle at a frequency non-synchronously related to the spin speed. The signal 
demodulated at this frequence is a measure of the scale factor unbalance of the summed pairs. The 
scale factor of an additional accelerometer is adjusted until a balance is obtained as illustrated by block 
diagram, Figure 1 5. 

Whether these scale factor balance loops will be required for the orbital gravity gradiometer depends 
on the amount of linear acceleration at the spin speed and twice spin speed angular modulation experienced 
aboard the spinning vehicle. Mechanisms for both type of distrubances are known and should their mag- 
nitude require it, these balance loops can be easily incorporated in the orbital gravity gradiometer. 

The status of the airborne gravity gradiometer development program is presented in Section 5 .0 of 
this report and the results form the foundation for the orbital gravity gradiometer. 
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Figure 14. Automatic Scale Factor Balance Loop No. 1 











Figure 15. Automatic Scale Factor Balance Loop No. 2 










4.0 ORBITAL VERSION OF ROTATING ACCELEROMETER GRAVITY GRADIOMETER 

The orbital version of the rotating accelerometer gravity gradiometer is discussed in this section. 
Miniature Electrostatic Accelerometers (MESA) are substituted for the Model VII used for the airborne 
gravity gradiometer and the spinning orbital vehicle serves as the rotating platform. The MESA’s potentially 
generate less noise to permit gravity gradient measurements to 0.03 EU in a ten second time period. The 
slov/ly spinning orbital vehicle (0.2 rad/sec) potentially provides an almost perfect rotating fixture as con- 
trasted with the bearing supported, torque motor driven device which is presently used and hence will 
eliminate many of the noise sources introduced by these components. 

The analysis covers the performance to be expected from a slowly spinning orbital gravity gradio- 
meter using four MESA’s. Two general areas of error mechanism are discussed, self generated noise from 
the MESA instruments and noise introduced by the rotating satellite. Critical requirements for the space- 
craft are outlined to allow gravity gradient measurements to 0.03 EU. 

4.1 MINIATURE ELECTROSTATIC ACCELEROMETER (FIGURE 16) 

The Bell MESA is of the force feedback type specifically designed for low level acceleration measure- 
ments in space applications. The MESA is a space qualified instrument and has been employed in a number 
of successful space missions including NASA’s Atmosphere Explorer and SERT 2 and the Air Force’s 573- 
5, Cannon Ball 1 and 2, and SPADES. 

Bell Aerospace’s experience in the operation of a number of instruments in which a capacitive type 
pickoff is employed, is that the instrument resolution is limited by self generated noise within the sensor 
itself and in the associated electronics. It is the self generated noise that is of prime concern in the selection 
of an accelerometer for a satellite version of a rotating accelerometer, gradiometer. It is not surprising that 
analysis shows that pendulous type accelerometers used successfully in a terrestrial rotating accelerometer 
gradiometer is not suited to the satellite version where 0.03 EU (10 sec time constant) or better perfor- 
mance is required. Preliminary analysis of noise from MESA’s indicate that instrument and system noise 
consistent with performance requirements can be met. 

Instrument and system self generated noise is analysed in Section 4.2 of this report, and indicate 
one, perhaps two, necessary modifications of the standard production unit. Further detailed analysis . 
might reveal some minor desirable dimensional changes of such parameters as gap lengths. None of these 
modifications impact on the basic qualified instrument design, a description of which follows. 

The proofmass, consisting of a thin-walled cylinder with a thin central flange, is electrostatically 
suspended.'by the application of voltages to the electrodes mounted to a carrier' fitted inside the cylinder. 
Suspended against any cross acceleration forces, the proofmass is free to move along the cylinder axis (the 
input or sensitive axis). 

Forcer assemblies are positioned on either side of the proofmass central flange. Each forcer assembly 
contains -three concentric ring electrodes. The central ring electrode is used in a capacitive pickoff- system 
to detect the position of the proofmass along the input axis. The two outer electrode rings are used to 
generate the rebalance forces to oppose the inertially induced forces from the accelerations to be measured. 

The operation of the force rebalance system is shown in block diagram form in Figure 1 7. 
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Figure 1 6. Bell Miniature Electrostatic Accelerometer (MESA) 
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The capacitive pickoff is energized from an oscillator typically operating at 100 kHz. The pick- 
off output signal which defines the proofmass position by amplitude and phase is amplified and demodu- 
lated. The signal is then suitably shaped for loop stability and the desired frequency response before 
being amplified for application to the forcer electrodes. The action of the loop is thus to constrain the 
proofmass along the input axis to the null of the pickoff system. 

Figure 1 8 (a) shows the most commonly implemented method of operating the capacitance pick- 
off system. The pickoff electrode to proofmass flange capacitances C 2 are connected in a bridge with fixed 
capacitances Q . The bridge is excited by a voltage Vp at 100 kHz. The bridge output is fed to a differential 
amplifier prior to demodulation. 

Figure 1 8 (b) shows the voltage applied to the forcer electrodes to generate the rebalance force. 

The reference voltage is continuously applied and is basically used to linearize the instrument output. 

The voltage V is the voltage applied from the forcer amplifier to satisfy the constrainraent loop. This 
voltage, or a multiple of the voltage, is the instrument output proportional to input acceleration. The 
voltage polarities are organized to maintain a net zero summed voltage with respect to the proofmass. 

Both Ej^ and V are dc voltages so that the instrument output is an analog dc voltage. 

The electrostatic suspension of the proofmass against across axis acceleration forces is usually ac- 
complished with a simple passive circuit. This suspension system can be scaled to permit instrument test 
on Earth in a Ig cross field. In a recent program, analog constrainment with a loop essentialy indentical 
to the input axis for the orthogonal cross axis has 'been developed to yield a 3 axis measuring instrument. 

For reasons explained in the noise analysis, analog constrainment in the cross axes will be necessary though 
not in a measurement mode. 

It is estimated that a single MESA complete with its associated electronics will occupy about 100 
in.^ weigh about 4 pounds and consume about 6 watts. 

4.2 Self Generated Instrument and System N oise 

This section is devoted to self generated noise, and does not consider noise and errors generated by 
other spacecraft running equipment or undesired spacecraft motions, which are covered in Section 4.3. 
However, one aspect of the environment, namely steady deceleration (drag) or acceleration of the spacecraft, 
does impact directly on the self generated noise and has to be considered here. Sources of self generated 
noise can be broken down into three categories. 

(1) True thermal noise (Brownian motion) in the MESA sensor unit. 

(2) Electronic noise in the circuits directly associated with the operation of the MESA. 

(3) Electronic noise in the circuitry associated with signal detection and interfacing with the 
spacecraft data transmission system. 

The preliminary analysis results given are for a 4 MESA unit system, with each sensor at a nominal 
radius of 30 cm. This yields a scaling sensitivity of 4.2 x 10"'^ g rms per EU for the cross gradient channel. 
The EU noise figures calculated are for a 10 sec time constant filter and apply to the lower sensitivity in- 
line gradient channel. The figures should be halved for the cross gradient channel. 
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Figure 18 . Pickoff and Forcing Arrangement 
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4.2.1 MESA Sensor Unit Thermal Noise 


Since the MESA is essentially a force measuring device designed to measure inertiaily induced ac- 
celeration forces acting of the proofmass, the instrument will indicate all forces acting on the proofmass 
including those not inertiaily induced. There are two mechanisms for thermally agitated forces acting 
directly on the proofmass within the sensor unit itself. These are: 

(1) Random collision of gas molecules with the proofmass. 

(2) Random electron flow in the proofmass reacting with any magnetic field which is present. 

The magnitude of these forces is related to the viscous damping coefficient associated with the 
mechanism; gas damping, and eddy current damping respectively for (1) and (2) above, by the following 
equation : 

fj.j^g = sj 4 kTc A f dyne (22) 

where 

k = Boltzmann’s constant in ergs/°K 

T = Temperature in °K 

c = Viscous damping coefficient in dyn/cm per sec 

A f = Bandwidth in Hz. 

The validity of this equation has been proven experimentally in testing conducted on pendulous 
type accelerometers used in the terrestrial version of the Bell Rotating Accelerometer Gradiometer. 

Equation (22) defines the total force noise, so that at a specific phase corresponding to one of the 
gradiometer output channels, and for a MESA proofmass weight of 1 gram, the acceleration noise from 
one sensor per channel is given by: 

1 

^ V^kTcAf 

^rms “ (23) 

981 

Since noise from four sensors adds on a rss basis, the noise from a pair of sensors increases to: 


®rms 


2v"4kTc Af _ y^SkTc Af 
98Iv^ 981 


(24) 


Ideally, the noise contribution from this source should be designed to be small with respect to 0.03 
EU. Solution of equation (24) for a noise contribution of 0.01 EU yields a maximum allowable value of 
c to be 3 X ICT® dyn/cm per sec, when the following parameters are used; 
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k = 1.38 xlCr^® 

T = 300 

A f = 0.1 rad/sec = 0.0 1 6 Hz ( 1 0 sec time constant filter) 

Srms“ ^ 

Very rough preliminary calculations indicate that to achieve this level of viscous damping from the 
gas contribution, the sensor unit internal gas pressure, assuming air, should not exceed 4x10"® torr. 

For any eddy current damping, the presence of a magnetic field oriented orthogonal to the input 
axis combined with a gradient along the input axis is required. There will be no eddy current damping 
from a perfectly uniform field so that it is the field gradient which is significant. Very rough preliminary 
calculations indicate that to achieve the desired level of eddy current damping, the field gradient should 
not exceed 0.8 gauss/cm. 

Restricting the magnetic field gradient to about 1 gauss/cm at the sensor location should not pre- 
sent any problems. Any problem that should arise in this area can be easily dealt with by adding magnetic 
shielding to the sensor or possibly by changing the proofmass material from beryllium to a higher electrical 
resistivity material. 

Achieving a sensor internal pressure of 4 x 10"® torr does raise a potential problem. It is believed 
that pressure this low could not be achieved and maintained in the standard hermetically sealed design. 

Two possible solutions that require further investigation are either venting the sensor in space or employing 
ion pumps in conjunction with sealed sensor units. 

4.2.2 Constrainment Loop Electronic Noise 


Electronic noise in the components of the constrainment loop causes forced motions of the proof- 
mass by the action of the closed loop. If these forced motions could take place without any resultant 
forces acting on the proofmass, then noise in the constrainment loop electronics would not cause noise in the 
output. There are, however, two sources of forces acting on the proofmass from imposed motions. They are: 

(1) Forces from negative spring rates to proofmass displacement due to electrostatics. 

(2) Forces from the proofmass inertia. 


Negative spring rates to proofmass motion along the input axis stem from voltages applied to both 
pickoff and forcer electrodes. The equation defining the spring rate is: 


dF 

dx 


1.76xlCr®AV^ 

-T dyn/cm 

d^ 


(25) 


where 


A = electrode area in cm^ 

V = applied voltage in volts rms 
d = gap length in cm. 
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Substituting the nominal values for a standard MESA yields tlie following: 

Pickoff spring rate = -5 x 10”^ g/cm per volt^ (26) 

Forcer spring rate = -1 x 1(T^ g/cm per volt^ (27) 

The inertia forces resulting from imposed motions can also be expressed as an effective spring rate 
at the frequency of interest. 


Inertia spring rate - - 


(2 n )2 

981 


g/cm 


(28) 


where 

S2 = rotation angular frequency in rad/sec. 

In order to assess the magnitude of the output noise from the constrainment loop electronic noise, 
it is necessary to establish the amplitude of the electronic noise referred to the preamplifier input. Knowing 
the pickoff gain at this point, the amplitude of the imposed motions can be calculated. It is then a matter 
of optimizing the operating voltages to minimize the output noise and determining any limiting environ- 
mental conditions. 

The most significant contributor to noise in the constrainment loop electronics is the first amplifier, 
usually referred to as the preamplifier. It is noise about the operating frequency of 100 kHz that is of con- 
cern in this stage. In subsequent stages it is the noise at the frequency of interest that is significant. Using 
the best low noise component currently available and assuming the preamplifier stage is operated with a 
gain of 80 volts/volt, the calculated rms voltage noise referred to the preamplifier input as a function of 
frequency as seen through a 0.016 Hz bandwidth is plotted in Figure 19. 

The gain through the pickoff bridge (at the input to the preamplifier) is known at 10 volts/ cm per 
volt of bridge excitation. 

The foregoing spring rate equations, constrainment electronic noise voltage ^d pickoff gain, permit 
calculation of output noise. There remains the optimum selection of two parameters, viz. the magnitude of 
the bridge excitation voltage Vp and the forcer reference voltage Ej^. The factors which have to be con- 
sidered in optimizing this selection are listed as follows: 


( 1 ) The pickoff bridge gain is proportional to— ^ while the pickoff electrode negative spring rate 


is proportional to — Therefore the ratio of 


spring rate 


is proportional to ——and for a 


given d, noise associated with the pickoff electrode spring rate is minimized by the selection 

of a low Vp. However, a low Vp means large motion since the motion is simply and 

^ gam 

these large motions are undesirable because of the forcer and inertia spring rates. 
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(2) The linear component of the rebalance force is proportional to Ej^V. Therefore the lower 
Ej^, the higher is V for a given g input. The voltage is continuously present resulting in 

a constant spring rate and constant magnitude noise contribution. The value of V, being 
proportional to the input g, generates a spring rate, proportional to and consequently a 
magnitude of noise proportional to the square of the input g. 

kg 

(3) The second order non-linearity error coefficient (in g/g^) is proportional to— — where 
kg is an asymmetry factor. 

(4) The dc bias Kq (in g) is proportional to Kg(Ej^)^ . Although the dc bias is not directly of concern 
in the gradiometer application, a parameter such as the bias temperature coefficient is also pro- 
portional to Ka (Er)^ and hence sensitivity at 2 to any 2 temperature variations. 

(5) The effective negative spring rate from the proofmass inertia is negative and therefore additive 
to the electrostatic spring rates. As the frequency of rotation is increased, this inertia spring 
rate becomes progressively more significant with respect to the electrostatic spring rates, re- 
quiring larger Vp for small motions for optimum conditions. 

Although the second order non-linearity and dc bias considerations cannot be ignored, preliminary 
investigation indicates that these aspects are unlikely to influence the selection of Ej^. It is noise that is 

by far the most significant factor influencing the selection of Ej^. 

A parametric study can be conducted to show the noise in the gradiometer output from electronic 
noise in the constrainment loop as a function of angular rotation frequency ^2 and g input (for example 
satellite drag and centripetal acceleration from rotation coupling in from tangential misalignment) for dif- 
ferent selections of Vp and Er- Two examples are shown in Figures 20 and 21, where the la EU plotted is 

the noise level as seen through a 10 second time constant filter. It is clear that from constrainment loop 
noise considerations, the lower the input g and the-rotation frequency, the lower will be the noise contri- 
bution. There is a secondary advantage here in that low input g and rotation rates permit the selection of 
low values for Vp and Ej^. Not only do these voltages generate a negative spring rate but also a net force 

if an asymmetry condition exists. Noise on Vp and Ej^ from the electronics generating these reference levels 
can potentially lead to problems if the nominals and asymmetry conditions are high. 

4.2.3 Suspension Electronics Noise 

The requirement to vritually eliminate gas damping in the input axis from thermal noise considera- 
tions, also eliminates damping in the two orthogonal suspension axes. A passive suspension scheme is 
therefore not practical. The proofmass will be suspended by the use of an analog constrainment loop 
similar in principle to the input axis constrainment loop. In this case a common electrode will be used for 
both the pickoff and forcing function. A MESA development program is currently underway to imple- 
ment such a scheme for the purpose of achieving a 3 axis measurement instrument. It has already been 
demonstrated that constrainment by this technique is practical. 

'The areas of concern for electronic noise in the suspension electronics are twofold: 

(1) Electronic noise in the suspension constrainment loop electronics will cause proofmass motions 
along the suspension axis. The potential exists for these motions to cross couple into the input 
axis pickoff system causing input axis motions. 
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Figure 21. Constrainment Loop Noise Vp = 0.25, Er = 0.05 
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(2) The application of pickoff and forcing voltages to the suspension electrodes potentially cause 
forces to cross couple into the input axis througli some asymmetry condition. In addition, the 
voltage application could generate a small negative spring rate along the input axis due to cer- 
tain geometric asymmetries. 

With regard to the motions along the suspension axis, since the pickoff gain is actually larger due to 
a shorter gap length, the motions should certainly be limited to no greater than those in the input axis due 
to its own constrainment loop. Even if as large a coupling as 1% existed (0.01 to 0.001% is more likely) this 
effect can be expected to be negligible. 

The only significant acceleration input against which the proofmass has to be suspended is the cen- 
tripetal acceleration from rotation. At of 0.1 rad/sec, this acceleration is roughly 3 x 1(T^ -g at a radius 
of 30 cm. Forcing voltages for this g input level will only be in the 1 or 2 volt region so that a negligible 
contribution to input axis stiffness will result. This assessment is based on experimental measurement with 
100 to 200 volt application required to suspend against the full earth’s g. The force noise in the suspension 
axis (10 second time constant) for a 3 x ICT”^ g input can be calculated by the same method used for the 
input axis, and can be shown to be roughly 2 x ICT^^ g. Cross coupling of forces from suspension axis to 
input axis is typically 3x10*® g/g so that the cross coupled noise is roughly 6 x ICT^® g which is negligible. 
Even at an = 1.0 rad/sec this cross coupled force noise increases by a factor of 100 to about 6 x ICT^. 
g per sensor or roughly 1 x ICr^^ g for a pair of sensors contributing 0.005 EU to the gradiometer noise. 

This preliminary analysis indicates that for angular rotation frequencies of 1 .0 rad/sec or less the 
contribution of noise from the suspension electronics is negligible. 

4.2,4 Signal Detection Electronic Noise 

The signai detection and processing electronics consists of the basic elements are shown in Figure 
13. The most significant noise contributor is the bandpass amplifier, this stage being at the lowest signal level 
point. High gain in this stage renders noise in subsequent stages negligible. Using the best low noise ampli- 
fiers currently available, the rms voltage noise at this amplifier input as seen through 10 sec time constant 
filter as a function of frequency is shown in Figure 22. 

The voltage scaling of the standard MESA output is a function of; 

( 1 ) Value of selected , 

(2) Max voltage swing capability of the constrainment loop forcer amplifier, the signal level 
attenuation before application to forcer electrode, and the maximum constrainment g level 
required. 

The noise contribution from the detection electronics as seen through a 10 second time constant 
filter as a function of rotation angular frequency is shown in Figure 23 for the following assumptions. 

Forcer Amplifier Swing 30 volts max. 

Signal level attenuation of 10 leading to a maximum constrainment capability of 6 x 10*^ g. 
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Since high, gain in the detection electronics renders noise in subsequent stages relatively insignificant, 
it is not expected that interface electronics will present a noise problem. 

4.3 GRADIOMETER SPACECRAFT INDUCED ERROR SOURCES 

For a perfect spinning system in a perfect environment tlxe only source of gradiometer error will be 
self generated noise. Compared to a terrestrial spinning mechanism in a Ig environment, the spinning orbital 
vehicle approaches the perfect system. However, there wUl be imperfections and some of these imperfections 
result in gravity gradient measurement errors. 

Without detailed knowledge of the spacecraft characteristics, a full error analysis is clearly not prac- 
tical at this time. Such an error analysis could only be conducted in a study phase in close cooperation with 
NASA and spacecraft personnel. However, there are a number of considerations relating to the spacecraft 
which are clear from the analysis of gradiometer errors. Some of these considerations may be specific to 
the Bell gradiometer employing rotating accelerometers, while some will apply to any gradiometer system. 
This section is devoted to a discussion of these considerations. 


In the consideration given to error sources, it should be borne in mind that the prime objective of 
the gradiometer is to measure the higher harmonic degrees of the earth and particularly haromincs in the 
band of about the 10th to 100th. Ideally, the gradiometer accuracy should be extended to the lowest 
possible harmonic to give a maximum overlap with harmonic data obtained from doppler tracking for 
correlation. From white noise considerations, the doppler tracking method is superior to the gradiometer 
at the lower harmonics of the earth and certainly below the 10th. The errors sources considered in this 
section are not white, and in general can be thought of as a bias output. An error mechanism generating 
a constant bias would be no concern. It is variations over an orbit of a bias producing error mechanism 
that is of importance. With an orbit time of about 90 minutes the 10th harmonic has a period of about 
9 minutes and the 100th of about 54 seconds. Bias variations outside this frequency range do not impact 
on the prime mission of the gradiometer. This is an important consideration and will be emphasized on 
more than one occassion in the following subsections. 

4.3.1 Spacecraft D rag 

Drag deceleration potentially impacts on the gradiometer performance in a number of ways. As 
shown in the preliminary analysis of self generated system noise, steady inputs must be limited to a few 
micro g. It is possible that this aspect is the most significant consideration with respect to drag deceleration. 

The drag deceleration will be seen in each MESA modulated at the spin frequency (in). Second 
order error coefficients potentially give rise to errors at the .gradiometer detection frequency of twice spin 
frequency (2^). 


Just considering the second order non-linearity coefficient K 2 , an error term is generated given by 


2 K 2 cos 2 52t 

1.2x Kr^° 


EU • ‘ ' 


where a^ = drag deceleration in g 


S K 2 = sum of second order non-linearity coefficients in g/g^ . 
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It will be seen that the error term contains a .^2 so if is limited to about 3 x 1 O'® g by selection 

of a suitable orbital attitude, 2 Ka would have to be very large before a significant error developed. The 
second order non-linearity coefficient is proportional to where is an asymmetry factor. Prior experi- 

Er3 

ence indicates that unless Ej^ is selected at an unreasonably low value, this error source will be negligible 
for values of a^j in the few micro g range. 

Asymmetry in the spacecraft profile can react with tlie drag deceleration to generate a spin frequency 
harmonic component in the drag deceleration and in a modulation of the spin frequency. Of concern is the 
1^2 component of deceleration and 2S2 modulation about the spin axis. 

The 1S2 content of the drag deceleration is modulated at 1112 by rotation and appears at the MESAs 
as a 212 acceleration. This signal is rejected in the summed output of the MESAs within the scale factor 
balance of the pairs. The automatic scale factor balancing described in Section 3.0, block diagram Figure 
14 is specifically designed to render this particular error mechanism to-be negligibly small. 

Angular modulation about the spin axis causes tangential accelerations to appear directly in the 
MESA input axes, and centripetal accelerations to couple into the MESA input axes by the tangential 
misalignment of the sensors. In the terrestrial version of the gradiometer this particular error source is 
minimized by a combination of a set up on test procedure and an automatic sum of pairs scale factor balance 
loop. This involves the deliberate introduction of a modulation about the spin axis at a discrete frequency 
close to but not harmonically related to the spin frequency. Such a scheme may not be beyond tlie wit of 
man in the orbital spacecraft application, but is clearly undesirable and in any case may well not be necessary. 

For a modulation angle about the spin axis of d peak radians, the tangential and centripetal accelera- 
tion is given by 

0.R. (my 

981 ^ 

The gradiometer error signal is given by 

k.g.R (2ny 
981 X 6 X 10-“ 

where k = rejection factor 

R = radius arm of accelerometer, cm 
f2 = spacecraft rotation frequency 

If a rejection factor k of 0.0025 is assumed together with an of 0.2 rad/sec then a 0.03 EU error will be 
generated by a modulation angle of 1.5 x 10“’ radians. 

The assumed rejection factor of 0.0025 is based on the accuracy of the scale factor balance between 
the summed pairs and on the net tangential alignment of the four sensors. The accuracy of the scale factor 
balance is dependent not only on the instrument scale factor in volts/g but also on the separation distance 
between the pairs of accelerometer sensors. Achieving instrument scale factor match to 2.5 parts in 1000 
is easy, and controlling the separation distance to that involves positioning of the sensors a nominal 60 cm 
apart to an accuracy of 1.5 mm. More difficult to achieve is the tangential alignment since this requires 
knowledge of center of rotation of the spacecraft. 
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Although not related to drag, the topic of this section, it is perhaps an appropriate point in the 
discussion to raise another requirement for tangential alignment. The constant rotation of the spacecraft 
generates a dc centripetal acceleration of about 1 .2 x ICT^ g at a radius arm of 30 cm for a 0.2 rad/sec rotation 
frequency. This dc acceleration will couple into the input axis by tangential misalignment and is thus a 
potential contributor to self generated noise as explained in Section 4.2.2. To limit the acceleration to the 
desired maximum of about 3x10"® g requires that the individual accelerometers be tangentially aligned to 
better than 0.0025 radians. This 0.0025 radian figure is consistent with the assumed rejection Ic of 0.0025 
above. 

Achievement of a 0.0025 radian tangential alignment in a 30 cm radius requires among other things, 
a knowledge of the center of rotation of the spacecraft to better than 0.75 mm. Prediction of the center of 
spacecraft roation to this accuracy is likely to prove impossible requiring that this problem be solved by a 
technique which does not involve rotation center prediction . 

One possible practical solution is to incorporate in the MESA mount an electromechanical mechanism 
for changing the tangential alignment. Such a mechanism could either rotate the sensor or translate the sensor 
tangentially. The output of each MESA will give an accurate indication of its tangential alignment by the 
magnitude of the dc level. This dc level will contain the instrument bias, expected to be very small, in ad- 
dition to the dc from tangential misalignment. With the spacecraft in orbit each MESA output can be moni- 
tored in turn, and on command the position servoed for a null dc output. This servo loop would not be 
continually activated because of noise generation, but activated at the start of the gravity gradient measure- 
ment mode. Should the center of rotation shift beyond 0.75 mm during the measurement mission (from 
fuel usage for example), the positioning system could be periodically reactivated on command. 

4.3.2 Spacecraft Temperature 

It is believed that temperature control of the MESA is highly undesirable from power usage con- 
siderations, and should not be necessary. Consideration has therefore to be given to the temperature en- 
vironment of the spacecraft. Errors associated with temperature variations occurring specifically at a 
frequency of 2i2, and to a lesser extent by have to be assessed in addition to the impact of the ambient 
temperature. 

Temperature variations at harmonics of the spin speed can certainly be expected from the radiant . 
heat of the sun acting on a slowly spinning spacecraft. Temperature variations of the MESA occurring at 
2S2 will react with the MESA bias temperature coefficient to generate a direct 252 error output. The bias 
temperature coefficient of the MESA is a function of the electrode voltages and for the proposed operating 
levels is predicted to be 2 x 10"^® per °C.’ Assumming uncorrelated coefficients, the error is given by 

2x2x1 0"^® 

^ EU per deg C at 252 

6 X 1{T“ 

= 6.7 EU per deg at 252 

A 0.03 EU bias would result from a 0.005°C variation at 252. A modulation of any 252 temperature 
variation with periods in the 54 sec to 9 minute range of this magnitude would not seem very likely. 

A temperature variation at 152 potentially impacts on the performance in two ways. First the 152 
temperature variation reacts with the MESA bias temperature to generate a 152 signal not related to scale 
factor imbalance. Such a signal causes, the automatic 152-scaIe factor balance loops (Block Diagram 
Fig. 14) to operate with an offset error. This error however, is assessed to be negligibly small. 
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Secondly and- more importantly the temperature variation reacts with the MESA scale factor 
temperature coefficient to change the scale factor at a 1^2 rate. This modulation reacts with the drag 
deceleration to generate a 2Q. error. If, for example, the' scale factor is modulated at m sin £2 1, then an mstru- 
ment seeing the drag deceleration modulated at sin £2 1, will have an output given by 

sin £2t 

1 + jn sin fZt 

= a^j sin £2t + 'A m a^ cos 2 £2t 

The uncompensated MESA scale factor temperature coefficient is of the order of 200 ppm per deg 
C. Assuming correlated coefficients, the error term without compensation, is given by 

4 X 14 X 2 X 10^"^ a^ 

— — ^ EUperdegCat li2 

6 X 1CT“ 

= 6.7 X 10® EU per deg C at 1^2 

For a drag deceleration of 3 x lO"® g this error term becomes 20 EU per deg C at li2. Open loop 
temperature compensation can be implemented to operate with probably a 5% accuracy bearing in mind 
temperature gradients that will occur at the 1J2 frequency. This would reduce the sensitivity to about 1 
EU per deg C at 1£2 yielding a 0.03 EU bias-for 0.03 deg C temperature variation. 

Probably thermal insulation in the mounting of the MESA from the spacecraft structure would 
considerably attenuate 112 temperature variations at the MESA itself. Again remembering that it is 
amplitude modulation of this 112 temperature variation in the 54 sec to 9 min periods that is of concern, it 
would seem likely that errors from this source will be small, but should be studied further. 

The only impact on gradiometer performance from ambient temperature changes is to change 
the gradiometer scale factor i.e., the volts per EU. This potentially impacts on the performance of the 
in-line gradient channel Wj^j^ - because of the relatively large standing level of about 3900 EU. Gradio- 
meter scale factor ambient temperature sensitivity can be expected from two main sources viz the scale factor 
sensitivity of the MESAs themselves and the temperature sensitivity of the radius arm. If the radius arm is 
part of the spacecraft structure, which would seem the most practical mechanical implementation, a tem- 
perature coefficient of up to 25 ppm per deg C is likely. Temperature compensation- of the MESA together 
with correction for the spacecraft temperature would probably reduce the overall scale factor coefficient 
to a few ppm per deg C. Temperature gradients would most likely be the factor which limits the accuracy. 
With a nominal gradient of 3900 EU an 8 ppm scale factor change results in a 0.03 EU error. If the ambient 
temperature within the spacecraft is expected to be ±I 0 deg C then variations with periods of 54 sec to 9 
mins might well be small. However, it is assumed that some form of temperature control is employed to 
maintain the spacecraft at ±10 deg C. The characteristics of any control will have to be considered for 
frequency content. 

4.3.3 Spacecraft Spin Speed, Attitude, Precession Rates, and Nutation 

As explained in Section 2.1, the minimum spin speed of the spacecraft is dictated by the data acquisi- 
tion rate and the highest harmonic sought. To acquire the lOOth harmonic degree requires 200 data points 
per orbit and with data acquisition at once per revolution this puts the minimum spin speed at a little over 
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0.2 rad/sec or roughly 2 rpin. Somewhat higher spin speeds, should they be desirable for other spacecraft 
reasons should not be ruled out. The biggest disadvantage for the gradiometer is the increased centripetal 
acceleration term. However, the proposed MESA tangential alignment mechanism, if adopted, should be 
able to tolerate larger centripetal accelerations than from the minimum acceptable spin speed. In some 
other error mechanisms, increased speed is an advantage e.g., 1^2 and 2S2 temperature variations. 

The proposed self generated demodulation reference signal scheme eliminates the requirement 
for any additional pitch attitude information, makes the actual spin speed non-critical and tolerates 
slow changes in spin speed which can be expected from aerodynamic and magnetic despin torques. 

The spin axis will be nominally orthogonal to the orbital plane, so that attitude-errors in roU and 
yaw have to be considered. In addition, angular rates about the roll and yaw axes generate acceleration 
gradients which the gradiometer promptly measures as it should, so that angular rates are a direct source of 
error. 


angle 


For roll and yaw attitude errors, the gradient signals are attenuated by cosine of the double error 


In Line Gradient =1 


W - W ' 

^ XX ** zz 


cos 2<f) 


Cross Gradient = 
where (f> = roll misalignment angle 


and \p - yaw misalignment angle 

Since the nominal value of is zero, the yaw error is not highly critical, but the nominal 3900 
EU in the in line gradient renders roll misalignments very sensitive. For a 0.03 EU error cos 2tJ/ has to be 
known to an accuracy of about 8 x 10"® . This suggests limiting the roll angle error to 0.12° or if the roll 
error is allowed to build up knowing the error angle to high accuracy - about 0.006° for 1° roll error, 0.003° 
for a 2° roll error. Again it is modulation of the roll attitude error in the 54 sec to 9 min period which is 
of concern, but in view of the high sensitivity this aspect requires further study. 


The angular rate sensitivity is simply given by 1 EU per IC® (rad/sec)^. A 0,03 EU error will be 
generated by an angular rate of 5.5 x lO”® rad/sec or about 1 .l°/hr. Consideration will have to be given 
to the attitude control system in the light of this sensitivity. 


Spacecraft nutation should be a true coning motion. Coning motion does not generate errors from 
the angular rates associated with the motion, so that nutation is not expected to be a problem area. 

4.3.4 Spacecraft Magnetic Environment 

The preliminary analysis of self generated noise indicated a requirement to limit the field gradient 
to 1 gauss/cm. The MESA performance is very insensitive to magnetic fields. A very small bias sensitivity 
can be predicted for a field orthogonal to the input axis with a gradient in amplitude along the input axis 
(same as self generated noise). Errors from this source should be negligible, but in any case can be easily 
dealt with by the addition of magnetic shielding to the MESA, if necessary. 
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4.3.5 Spacecraft Vibration 


The potential exists for the MESAs to experience acceleration inputs from the vibrations transmitted 
through the spacecraft structure from other equipment aboard the spacecraft. Common mode inputs to all 
four MESA are rejected within scale factor balance and alignment accuracies. However, some equipment 
such as horizon scanners for example are likely to produce inputs which are not common mode. Since any 
vibration and noise, producing equipment will be fixed with respect to the accelerometers of primary concern 
is the magnitude of vibration at and near (about 0.4 rad/sec). 

A vibration power spectrum which is flat in the region of 2Q, which couples into one accelerometer 
ONLY, of magnitude per rad/sec will cause a gradiometer noise output of 

(4.2tl0-) 

Thus for a 0.01 EU^/rad/sec accuracy (0.03 EU in 10 seconds), the value of 'S' representing a non- 
common mode .input cannot exceed about 18 x 10"^ per rad/sec. Although this is a very small density 
amplitude it is .applicable to a very low frequency of 0.4 rad/sec. It is clear, however that further study in 
this area is required. -Operating in the gravity gradient measurement mode between the poles with all noise 
producing spacecraft equipment shutdown may or may not be a feasible approach. Mounting of the ac- 
celerometers to a dedicated structure with careful design in the mounting of this structure to the spacecraft 
may well substantially reduce non-common mode inputs. 

4.3.6 Spacecraft Orbital Height, Spin Speed and Attitude Corrections 

The effect of aerodynamics and magnetic fields will require that periodic thrusting be applied to 
maintain orbital height, spin speed and attitude. The foregoing consideration of error sources suggest that 
performance degredation might be experienced when thrusting is activated. Passes over the earth poles will 
be made in each orbit. Use may be made of this fact in two possible ways. One is to use data acquired in 
each pass in a calibration mode. The other is to disregard gradient data acquisition while required thrusting 
is .accomplished. Depending on how often such thrusting corrections are necessary, possibly a combination 
of the two advantages could be beneficially implemented. 
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5.0 STATUS OF ROTATING ACCELEROMETER GRAVITY GRADIOMETER 
DEVELOPMENT PROGRAM FOR SAMSO 


A status report on the feasibility demonstration model of the rotating accelerometer is added to 
this report to substantiate the soundness of the approach and lend a high confidence factor of success to 
the orbital gravity gradiometer. The objective of the airborne feasibility model is to demonstrate a per- 
formance of 1 EU as measured in all attitudes and under dynamic conditions. The attitude most closely 
approaching the orbital case is a slightly off vertical rotation axis (minutes of arc) causing a small amount 
of acceleration to couple into the sensitive axis of the accelerometer. 

For the initial test program a GFE precision rate table was used to furnish the rotating fixture. The 
photograph. Figure 24, illustrates the setup. The four Model VII accelerometers can be seen mormted to 
the instrument block. The constrainment electronics and the detection electronics is mounted on vector 
board with each component available for troubleshooting arid investigative testing 

In the background is the 100 kgm gradient inducer which can be traversed on an overhead rack for 
calibration purposes. 

The block diagram, Figure 25, illustrates the signal flow. The four accelerometers aj through 34 
are schematically shown with their constrainment electronics through on the rotating fixture 

indicated by the circle. The four -constrainment currents are summed and amplified by the detection 
electronics prior to being passed through slip rings. The output of the detection electronics is demodulated 
at sine 20 1 and cos 20 1 and passed through ten second filters. These are the outputs of the gradiometer 
which are displayed on an analog strip chart and digitally tape recorded for computer analysis. The output 
of the detection electronics is demodulated also at sin Ot, cos Ot and sin cjgt to establish the scale factor 

unbalance of the accelerometers. These signals are integrated and currents generated into the scale factor 
adjustment coils. (Og is the angular modulation frequency imposed in the spin speed of SI to permit detec- 
tion of scale factor unbalances between the summed pairs of accelerometers. 

The output of individual -accelerometers are also passed through slip rings for diagnostic purposes. 
One accelerometer output is demodulated at sin f2t and cos S 2 t to indicate tlte off vertical tilt of the 
rotation axis. 

The chart on Figure 26 gives the typical performance figures of an overnight run. An HP 830 
computer continuously records the output of both channels of the gradiometer and on-line analyses 128 
consecutive 10 second average data points for mean (bias, trend (change of bias with time)) and random- 
ness around a best fit straight line. 

The randomness of 80 blocks of 1 28 data points is 1 .07 EU with the worst reading and best reading 
not exceeding ± 20 % of that value. 

The bias of 40 consecutive data blocks exhibits a one sigma randomness of less than 0.4 EU with 
negligible trend for an overnight run. The' overnight bias stability for the 40 data blocks of each gradio- 
meter output channel is plotted on Figure 27. The excellent bias characteristics for an overnight run with- 
out temperature control is indicative of the soundness of the approach. 
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Figure 24. Rotating Accelerometer Gravity Gradiometer 
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ROTATING ACCELEROMETER GRAVITY GRADIOMETER 

ROTATION SPEED: 3/8 Hz 

OFF VERTICAL TILT OF ROTATION AXIS: APPROX 2 ARC MINUTES 

DATA ACQUIRED: THROUGH DOUBLE SECTION 10 SEC FILTER RECORDED APPROXIMATELY 

EVERY 10 SEC 

DATA ANALYZED: EVERY 128 DATA POINTS (APPROX 1280 SEC) 40 SEQUENTIAL BLOCKS OF DATA 

FOR BIAS, TREND AND ^a RANDOMNESS BOTH GRADIOMETER CHANNELS 

RESULTS: la RANDOMNESS: MEAN 80 DATA BLOCKS . 1.07 EU 

WORST 80 DATA BLOCKS 1.27 EU 

BEST 80 DATA BLOCKS 0.87 EU 

BIAS: SINE CHANNEL - MEAN 40 DATA BLOCKS 23.75 EU 

la DEV 40 DATA BLOCKS 0.39 EU 
TREND 40 DATA BLOCKS 2.2 x 10~3 EU/HR 

COS CHANNEL - MEAN 40 DATA BLOCKS 7.54 EU 
la DEV 40 DATA BLOCKS 0.25 EU 
TREND 40 DATA BLOCKS 1.57 x 10"3 EU/HR 

NOTE: APPROXIMATELY 100 EU BIAS BUCKED OUT BY PLACEMENT 100 KG M LEAD BALL 


OJ 


Figure 26. Performance Summary 
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The performance of the rotating accelerometer gravity gradiometer using Model VII accelerometers 
and the Holloman rate table is limited by the noise of the Model VIIs (non-rotating) v/hich has been meas- 
ured to be in the order of 0.9 EU and a small amount of additional noise introduced by the rate fixture. 

The effectiveness of the common mode rejection of external disturbances using four Model VII 
accelerometers in combination with the automatic scale factor loops is clearly illustrated by the noise 
measured along the various axes of the Holloman rate table expressed in EUs as illustrated on Figure 28. 
The anticipation of being able to measure gradients to an accuracy of 0.03 EU is based on the use of the 
lower noise MESAs and the much smoother rotation provided by a spinning orbital vehicle. 
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T a RANDOMNESS ABOUT SPIN AXIS AT f2 = 1/2 Hz : 11,800 EU 


Figure 28. Nominal Amplitudes of 2Q, Experienced by a Sensor 



6.0 CONCLUSIONS AND RECOMMENDATIONS 


6.1 CONCLUSIONS 

The following conclusions have been drawn based on this study. 

• Gravity gradients to an accuracy of 0.03 EU in ten second average are possible with a slowly 
spinning orbital vehicle using four modified versions of the space qualified Miniature Electro- 
static Accelerometer. 

• At an orbital altitude of 300 km and a spin rate of 2 rpm, global gravity mapping to the 80 
harmonic degree will be feasible with a 0.03 EU gravity gradient measuring accuracy. 

• The spinning orbital gravity gradiometer self generates the demodulation reference and hence 
all signal processing to obtain the gradiometer outputs can be conducted onboard the satellite. 
Transmission of one 8-bit data point for the gradiometer channels every 1 5 seconds suffices 
for global gravity mapping and places only the most moderate requirements on data trans- 
mittal and handling. It is planned to buck out the nominal vertical gravity gradient of 

4000 EU onboard the orbital vehicle. 

• The attitude of the spin axis of the satellite should be maintained to 0.12 degrees in roil with 
respect to vertical to the orbital plane. Alternately knowledge of the roll attitude angle is 
satisfactory. The yaw attitude accuracy of the spin axis is not critical. 

• Angular input rates about the yaw and roll axis should be under 1.10 deg/hr. Care will have to 
be exercised in the design of the space craft to avoid angular jitter induced by other equip- 
ment on the space craft from exceeding that rate. Attitude correction over the polar regions 
where much redundant data acquisition occurs is one option to be considered. 

6.2 RECOMMENDATIONS 

The two critical areas on which feasibility of an orbital gravity gradiometer depends are: 

• The noise of the modified MESAs - The noise of a pair of MESAs modified for this applica- 
tion is readily tested in the laboratory and it is recommended that such a test program be 
conducted. 

• The spinning orbital vehicle requirements are outlined in this report for the gravity gradio- 
meter mission. 

■ A review of a preliminary layout of the GRAVSAT with Fairchild and NASA Goddard personnel 
indicate it to be suitable for an orbital gravity gradient mission. A detail analysis with appropriate space 
craft engineering personnel is recommended to arrive at an optimum configuration for gravity gradient 
mission. An alternate possibility is a Scout Launcher vehicle. This should be explored further because 
of the potential cost effectiveness. 
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